EXISTENCE OF A DEGENERATE SINGULARITY IN 
THE HIGH ACTIVATION ENERGY LIMIT OF A 
REACTION-DIFFUSION EQUATION 
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Abstract. Wc consider the singular perturbation problem 

Qi^ ' where (3e{s) ~ 7/3(f)j /3 is a Lipschitz continuous function such 

<C ■ that /? > in (0, 1), /? = outside (0, 1) and (3{s) ds = i. 

. We construct an example exhibiting a degenerate singularity as 

, efc \ 0. More precisely, there is a sequence of solutions u^^ — » u 

as A; — > cx), and there exists x° G 9{u > 0} such that 

— ^ — as r ^ 0. 

r 

Known results suggest that this singularity must be unstable, which 
, makes it hard to capture analytically and numerically. Our re- 

sult answers a question raised by Jean-Michel Roquejoffre at the 
FBP'08 in Stockholm. 



1. Introduction 
In this paper we are concerned with the singular perturbation prob- 



H ' lem 



(1.1) Au, = p,{u,) in n, 

where f2 is a bounded domain in M", e > and Pe{s) = 7/5(f). Here 
/3 is a Lipschitz continuous function such that /3 > in (0, 1), /3 = 
outside (0, 1) and Jq f3{s) ds = ^. 
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This problem arises in the mathematical analysis of equidiffusional 
flames (see [3], [1]), in which case e is proportional to the inverse of the 
activation energy. 

Formally, as e ^ 0, the solutions converge to a solution u of the 
free boundary problem 

An = inl]\9{n>0}, 
^ ' ' u = 0, - {u-f = 1 on d{u > 0}. 

On a rigorous level, L. Caffarelli established a locally uniformly Lips- 
chitz estimate for bounded solutions {u^} (see [5]). In [10], C. Lederman 
and N. Wolanski proved that m is a viscosity solution of (11.21) . They 
also proved that u satisfies the free boundary condition in a pointwise 
sense at non-degenerate free boundary points at which there is an in- 
ward unit normal of {u > 0} in the measure theoretic sense. Here 
"non-degenerate" means that 

^-n-i / u"*" > c > for r < tq. 

JBrix) 

A related convergence result has been proved for the p-Laplace operator 
in [6]. 

From [H Lemma 3.4] we know that all free boundary points of local 
minimizers are non-degenerate. This suggests that if there are degen- 
erate free boundary points of a limit function u, then they must be 
unstable. 

It has been known that cross-shaped free boundaries occur for do- 
main variation solutions (see the Introduction of [12]), for example 

u{xi,X2) = \xl — xl\ (cf. Figure[I]). 

John Andersson suggested a triple junction example (see Figure [2]) 
where the solution is homogeneous of degree 3/2; domain variation 
solutions with a degree of homogeneity G (1, 3/2) are not possible. 
However it is not so obvious whether these examples arise as limits of 
the reaction-diffusion equation (11.11) . This is the question Jean- Michel 
Roquejoffre posed at the FBP'08 in Stockholm. The main result of the 
present paper gives a partial answer to that question: 
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Figure 1. Cross-like Singularity 




Figure 2. Triple Junction Singularity 

Theorem A. There exist uniformly bounded solutions u^^ : -Bi(O) C 
^ M such that 0, u^^ u ask ^ oo, and there is a degenerate 
free boundary point of u in i?i(0). More precisely, there is G -Bi(O) fl 
d{u > 0} such that 

u(x° + rx) ^ „ ^2 
lim = U for every x G M . 

r^O r 

Remark 16.51 shows that we may construct degenerate points with an 
arbitrarily high number of symmetry lines. See Figure [3] for an example 
with 4 symmetry lines. 

Section [3] also suggests a numerical approach for the construction of 
this unstable solution which would be otherwise hard to capture. 

2. Notation 

Throughout this article will be equipped with the Euclidean inner 
product X ■ y and Br{x'^) will denote the open ball of center x^ with 
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Figure 3. Imaginable e-solution 

radius r. When not specified, is assumed to be 0. We will use 
the polar coordinates (r, 6) in R^. We use the ra-dimensional Lebesgue 
measure and the m— dimensional Hausdorff measure Ti™. We shall 
often use abbreviations for inverse images like {u > 0} := {x E Q : 
u{x) > 0}. 



3. Construction of e-soLUTioNS 

We are going to use a method introduced in [2]. 

Let K = {{r cos 6', r sin 6) : < r < 1,0 < 6 < vr/m}, where m is an 
integer > 3. For g e C^^dBi n dK), ao > « > 0, C^iK) will denote 
the subspace of C°(A') consisting of all functions with boundary values 
g on dBi fl dK. We consider the following problem: 



An(x) = I3^{u{x) - u(0) + e) in K, 

(3.1) u = g ondKndBi, 

1^ = on dK\dBi. 

Let T = T,^g : C^{K) ^^"(i?) be the operator defined by 

AT(m) = p,{u{x) - u{0) + e) in K, 

T{u)=g ondKndBi, 

^ = ondK\ dBi. 

Using Schauder's fixed point theorem, we can prove the existence of 
a solution of equation (13.11) : 
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Proposition 3.1. For each g G C°'°{dBi fl dK), T has a fixed point 
Ue G Cg{K) for some a > depending only on ao. 

Proof. Let v G C^(i^), F{x) = /3,{v{x) - v{0) + e). Then there is a 
iy^'^(-ft')-solution u of the problem: 

Au = F in K, 

u = g on dK fl dBi , 

1^ = ondK\dBi. 

We extend m to a function u : Bi \ {0} ^ M by even reflection: 



u{r, Tr/m + 6) = u{r, vr/m — 0), if < < n/m, 

u{r, jir/m + 6) = u{r, jir/m — 6), if < 6 < ir/m, 

for j = 1, 2, ... , 2m — 1. It follows that £t is a solution of 

Au = F mBi\ {0}, 
u = g on dBi, 

where F and g are defined by even reflection. As the origin is a set of 
vanishing 2-capacity (see [7]), m is a weak solution of Am = F in _Bi. 

Applying the regularity theory of elliptic equations (see for example 
P Corollary 9.29]), we see that u G C°'{Bi) for small a, and that 
I l^^l |co>"{Bi) ^ C, where a G (0, ao) and C are constants depending on 
Q^O; llfi'llcO'^o and | |-F| |ioo(^). Hence T is a continuous compact linear 
operator from Cg{K) into itself for small a > 0, and 

\\Te,g{w)\\c^^K) < C 

for every a G [0,1] and every solution of the equation 

Aw = (7(3,{w{x) -w{0) + e) in K, 

w = g on dK fl dBi, 

1^ = ondK\dBi; 

here C is a constant depending only on g, (3 and e. 

From Schauder's fixed point theorem (see [U Chapter 11]) we infer 
that T^^g has a fixed point G Cg{K). □ 
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4. Convergence 

From now on all of our functions are defined in -Bi(O) by the above 
reflection. We choose one non-constant g G C'^{B). 
Let fe(x) = — Me(0) + €. Then Ve is a solution of 

Av,(x) = (3,{v^{x)) in 5i, 
t;,(0) = e, 

Ve = g — Ue{0) + e on (95i. 

By the maximum principle, |Me(0)| < 1 |ioo(B^), therefore | |we| |2,oo(5^) < 
C, where C depends only on 1 |loo(5^). 

First we state a result from [5] proving a uniform Lipschitz estimate. 

Proposition 4.1. Let be a family of solutions to Aw^ = l3^{w^) in 
a domain Q G such that | |L°°(n) ^ C for some C < oo. Let 
K G Q be a compact set and let r > be such that Br{x^) C Q for 
every x° G K . Then there exists a constant L = L{t, C), such that 

\VWt,{x)\ < L ioT X E K. 

If the boundary data is smooth enough, we have the following gra- 
dient estimate near the boundary (see [9]). 

Proposition 4.2. Let Q C be a bounded domain with smooth 
boundary and let G C'^{Q) fl C^{Q) be a solution of 

Awe{x) = Peiwtix)) in Q, 
We = f on dfl, 

where f G (7^(0), ||/||c2(n) < Ci. Then there exist eo > 0, C < oo 
such that 

\Vw,{x)\ < C{1 + \\ogd{x,dn)\) 

for e G (0,eo), x G Q, where d{x,dQ) is the distance from x to the 
boundary of Q . 

Now we are ready to prove the convergence of {ve^}. 

Proposition 4.3. There exist a sequence {v^^} and a locally Lipschitz 
continuous function vq such that 
1) vq uniformly on Bi, 
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3) vq is harmonic in Bi \d{vQ > 0}, 

4) Av^- ^ fi as measures on Bi; here fi is a locally finite non-negative 
measure supported on the free boundary d{vQ > 0}. Therefore 

Avq = fi in Bi, 

that is 



Vvo -Vcp dx+ (p dn = 

Bi Jbi 

for every </> G Cq(-Bi). 

5) Vo{r, jn /m + 6) = Vo{r,jTT/m — 6) forO < 6 < ir/m and 1 < j < 2m. 

Proof. 1) By Proposition 14.21 {v^} is equicontinuous. Since {v^} is uni- 
formly bounded in C°(i?i), by the Arzela-Ascoli theorem there is a 
sequence {v^.} that converges to some vq in C^{Bi). By Proposition 
14.11 it is easy to see that vq is locally Lipschitz continuous. 

2) For < ?7 G C^{Bi), we obtain from the local Lipschitz bound 
that 

(4.1) / ?7|VfoP dx < liminf / ?7|Vfe ^ dx. 

Jbi ^'^0 Jbi 

If we prove the converse inequality, the result will follow. Multiplying 
the ej-equation (II. ip by rfv^., we obtain that 



v^yrj ■ Vw,, dx+ vl^^e.l dx = ~ r]v^^/3^X^^J dx < 0. 

Bi Jbi Jbi 

Letting e-j ^0, we get 

(4.2) limsup / ?7|Vfe-|^ dx < — vqS/t] ■ Wvq dx. 

ei->0 Jbi Jbi 

For each 5 > 0, we define 

vo{x) — 5, if vq{x) > (5, 

vsix) = { vo{x) + 5, if vo{x) < -5, 

otherwise. 
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Multiplying the ej-equation by rivs and integrating, we have 

vsVrj ■ Vt'ei dx + rjVvs ■ Vw^. dx = 0. 
Bi Jbi 

Letting ej — > 0, it follows that 

vsVr] ■ Vvo dx + r]\'Vvs\'^ dx = 0, 
Bi Jbi 



and letting 5 — > 0, we get 



VoVt] ■ Vwo dx + '7|Vfo| dx = 0. 
Bi Jbi 

Using (14.20 . we obtain 

(4.3) limsup / ?7|Vfe. p dx < ?7|VfoP dx. 

ei-*0 Jbi J Bi 

Combining (14.11) and (14. 3p . we obtain 

(4.4) lim / ri\Vve-\^dx= / ri\SJvQ\^ dx., 
^'^0 Jb^ ' JSi 

which implies that 

7]^''^Wv,^ 7]^^^Vvo in L\Bi). 

For 7] such that t] = 1 in D GG Q, it follows that 

Vw,, Vvo in L2(L'). 

Therefore — > fo strongly in ly^'^(D). 

By a diagonal sequence argument, assertion 2) of the proposition 
follows. 

3) Since vq is continuous, the sets {vq > 0} and {vq < 0} are open. 
Let x'^ G {vq > 0}. From the fact that v^. vq uniformly on Bi, there 
exist r > and N eN such that ^^.(x) > vo{x^)/2 > for x G Br{x^), 
i > N. Thus v^- is harmonic in Br{x^) for i > N and the same fact 
holds for t'o- In the same way we prove that vq is superharmonic in 
{^'o ^ 0}°. On the other hand, f^. is subharmonic in Bi, so that vo is 
also subharmonic. Therefore Vq is harmonic in {vq < 0}° and 3) follows. 
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4) Multiplying equation A?;^^ = P^.{veJ by G C^{Bi) and inte- 
grating by parts, we get 

/ Vt;,, ■ V0 rfx + / (3eXvei)(p dx = 0. 

Jbi Jbi 

We know that f^. — > vq in W^i^^{Bi), hence 

(4.5) / Vv,^ -Vcj) dx^ I Vvq ■ V(f) dx. 

Jbi Jbi 

On the other hand, 

[ PUveJ^dx<C. 

JBi 

This L^-bound implies that there exists a locally finite non-negative 
measure fi such that for a subsequence which we still call {fe^}, Pti{v^.) — ; 
fi as measures in Bi. Passing to the limit in fl4.5p . we get 



Wvq ■ V0 dx + (p dfx = 0, 
Bi Jbi 

which implies 

At>o = fi in 

Since we know that vq is harmonic in Bi \ d{vQ > 0}, we conclude that 
supp fj, C d{vo > 0}. 

5) follows from the fact that v^. has been defined via even refiection. 
This property is preserved by the uniform limit vq. □ 

Let Bt,{z) = Jq Pe{s) ds and Xe{x) = 2B^{v^{x)). It follows that 
^ ^ Xe{x) < 1, and the uniformly Lipschitz estimate for implies the 
relative compactness of Xe in Aoc(-^i)- 

Proposition 4.4. {Xe} is precompact in L^{D) for each D CC Bi. 

Proof. By Proposition l4.lt there exists an L > such that 1 1 Vf^l |loo(£)) < 
L, hence 

/ |VXe| dx = 2 \/3,{v,)\/v,\ dx<2L / (3,{v,) dx < C. 
Jd Jd Jd 

Therefore {xe} is bounded in W^'^{D). By the Sobolev embedding 

theorem {Xe} is precompact in L^{D). □ 
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Thus we may assume that Xeii^) Xo{^) locally in L^. Similar to 
T2I Lemma 4.1], we have 



Proposition 4.5. Xo{^) ^ {0? 1} foi" o-e. x E B 



1- 



Proof. For 5 G (0, 1/4), B{s) = 2 p{t) dt and As = {x e D CC Bi : 
Xo G [25, 1 — 25]} we have 

C\As) < C\As n {xe^ ^[S,l- S]}) + C\As n {S(^) e [5, 1 - S]}) 

< c\D n - xol > 5}) + n G [^-^(5), ^-^(i - 5)]}) 

< n {\x.^ - xol > 5}) + Ci(5,/5)e„ / dx 

JD 

< C\D n - Xol > 5}) + C2(5, /3)e™ ^ 

as m ^ cxD. □ 

5. MONOTONICITY FORMULA 

In this section we state a monotonicity formula proved in [13] , which 
is a key tool in constructing a degenerate singularity. 

Proposition 5.1. Let x° G -Bi, and 

Then $e satisfies the monotonicity formula 

$e(cr)-$,(p) > / /■ 2{Vv,-u-^^^f dV}dr forO<p<a. 



Letting e ^ 0, 



^[r)=r-^ \Vvo\'^ + r~^ Xo - Vq^ dV} 

JBrix") JBrixO) JdBrixO) 

satisfies the monotonicity formula 

$((t)-$(p) > r r^^ [ 2{Vvo-u-^^y dV}dr forO<p<a. 
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6. Existence of a Degenerate Free Boundary Point 

In this section we are going to prove the existence of a degenerate 
singular free boundary point for v^. We start with some lemmas. 

Lemma 6.1. Let w be a solution of 

Aw = /3,{w) znncW and let Xe{x) = 2B,{w{x)). 

Then 

|Vty|^div — 2VwD(j)Vw + Xediv (p dx = 



'n 



for (j) e C^,{n; 
Proof. Integrating by parts, we get 



I Vw| div — 2S/wD(j)SJw + Xediv dx 

n 



^ dxj ^ dxi dxi dxj 

2 [ "^^^ ■ ^ f ^ + /9 ( ) V (j) d 

JqA^ ^ dxi dxidxj ^ ^ dxi \dxidxjj 



' i,j=l ■' i,j=l 

2 Vw- (j){Aw - (3e{w)) dx 



n 



0. 



Lemma 6.2. If Vq = in Br for some r > 0, then xo = ^ in- Br- 
Proof. By Lemma [6.11 we have 

/ \Vv,fdiv - 2Vv,^D(jNv^^ + Xe.div = 



□ 
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Figure 4. The set Pi 

for (p G C^{Br] M^). Letting ej — *■ 0, we obtain from the strong conver- 
gence of that 

= / I Vfopdiv — 2Vfo-D0Vfo + Xodiv 
Xodiv dx, 

hence Xo = const in B,.. We know that Xo{^) ^ {0, 1} for a.e. x G 
therefore Xoi'-^) = or Xo{^) = 1 in i^r- 
For < c < 1/4, let Pj = {x e 5r : cej < v^X^) < (1 - c)ei} and 
= {x ^ Br : Vei{x) > cei}. Supposing towards a contradiction that 
Xo = in Br, and recalhng that Xei Xo in L^{Br) we have 

C'iQ.) < C'{{x : Xe. > 2B{c)}) ^ 

as z — > oo. 

Let 5 G (0,r/4), 4 = {t : 5 < t < r, min > 6^/2} and 

4^. = [5,r] \ 4. Then we have C\l,^) < C'^{Qi)/6 ^ as i ^ 00, 
hence there exists an eo > such that > r/2 for < eo- As 

v^- is subharmonic in Bi, ej = "^£^(0) < maxv^- for every t G [0,r). It 

dBt 

follows that for each t G I^., there exists a point G dBt such that 
= ej/2 (cf. Figure H]). By the uniformly Lipschitz continuity 
of v^^ there exists a. d > independent of e, such that Bdei{x\) C Pj. 
Therefore 
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On the other hand, since Awe^ — >^ in Br, 

C\Pd/e.<C [ PM^O. 



JBr 

Thus we get a contradiction and therefore xo = 1 in -Br- D 

Lemma 6.3. Let Vo{x^) = for some x° G Bi, OKr^^Oask^oo 
and VokM = !^oi2-±!M) _ Then there exists a blow-up limit Vqq such that 
for a subsequence the following holds: 

1) VQk vqq in Ci°'"(M2) for every < a < 1, Vwofc ^ Vt^oo weakly* 
tnir^M'), 

3) vqq is homogeneous of degree 1. 

Proof. 1) Since Vvok{x) = Vvq^Xq + r^x), {vq^} is locally uniformly 
Lipschitz continuous in and 1) follows. 

2) Let D CC M^. Then there exists a C < cxo such that |Vfo/c(2;)| < 
C in D, hence Vok is bounded W^''^{D). It follows that there is a 
subsequence, which we still call vok, such that vok vqq weakly in 
W^'^{D) and 

||'?^oo||vFi'2(z)) < liminf ||fofc||i4/i,2(B). 

k — ^oo 

Since A^o = and supp ji C d{vQ > 0}, we get 

i;ofc(x)Awofc(x) = 0. 
Let T] e Co°°(R2). Then /j^^ r/| Vt^ofeP = - ^2 ^^ofcV?/ ■ VvQk and 

limsup / ri\S/vok\^<- I vooWrj ■ Vvqq. 

For each 5 > 0, we define 

Voo{x) - 5, if vqq{x) > 6, 
vs{x) = <; voo{x) + 5, if voo{x) < -S, 
otherwise. 

Multiplying the equation for fofc by r]Vs, we get 

ri{x)vsix)Avokix) = 0, 



/ 
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if k is large enough. Hence 

r]Vvs ■ Vvok = - v&S/ri ■ Vvok- 
Letting — oo, it follows that 

rjVvs ■ Vvoo = - vsVr] ■ Vvoo, 
and letting 5 ^ 0, we get 

/ '7|Vfoor = -/ fooV?7-Vfoo- 

Hence 

limsup / r7|Vt;ofcP < / 77|Vt'ooP 
for every 7] G ^^(M^). It follows that Vok -> woo in W^'^{D). 

3) For < i? < S < oo, 

J Br JOBr 

Since vq{x^) = and vq is Lipschitz continuous, $(r) is bounded. Con- 
sequently we obtain from the monotonicity formula Proposition 15.11 
that 

2|x|"^(Vfofc(x) ■ X - vok{x)fdx. 



'Bs\Br 

Letting /c oo, we obtain that Vvoo{x) ■ x = Vqq^x) a.e. in M^, hence 
Vqo is homogeneous of degree 1. □ 

Now we are ready to prove Theorem A of the Introduction. 

Theorem 6.4. There is a free boundary point x^ G d{vQ > 0} such 
that Vq is degenerate at x^. More precisely, 

VQ(x° + rx) 
lim = 

r— ►0 r 

for every x G M^. 
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Figure 5. Touching the Free Boundary 



Proof. Case 1:0^ d{vo > 0}. 

There exists an r > such that Vq{x) < in B^. Since vq{0) = 0, 
the subharmonicity of vq imphes that vq = in Br- Suppose that r is 
the largest number such that vo = in Br. Noticing that r < 1 as the 
boundary values of vq are not constant. It follows that there is a point 
x° G dBr such that x° G d{vQ > 0} (cf. Figure E]). 

By a translation and a rotation we may assume that x° = and 
Wo = in Br{—r,0). Let vqq be a blow-up limit of vq at i.e. the 
limit of ""^^rf'''^^ as — s> oo. We are going to show that vqq = 0. 

Suppose towards a contradiction that vqo ^ and let S = {vqq > 0}. 
Then vqo is harmonic in S. Moreover, we know from Lemma 16.31 that 
Vqo is homogeneous of degree 1, so that, solving the resulting ODE for 
Vqo, Vqo = cx ■ u in S for some c G M and z/ G dBi. Noticing that 
Vqq = in the left half plane, we conclude that Vqq = cmax(0,x ■ u), 
where c > and u = (1, 0). 

Let Xr{x) = Xo(^a^) and let xoo be the limit of Xr-k of the above 
sequence k oo. By Lemma [6l2l xo = 1 in r, 0), hence xoo = 1 
in the left half plane. Moreover we know that Xo{^) = 1 if vo{x) > 0, 
hence xoo = 1 in the right half plane. Thus Xoo{x) = 1 in M^. 



From 








for every (/) G C^(M' 



M^), we infer that 
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for every G C^(M^; M^). This however contradicts foo = cmax(0,a; ■ 
u) . We conclude that the blow-up limit voo must be the constant func- 
tion in M^. 

Case 2: G d{vo > 0}. 

Suppose towards a contradiction that a blow-up limit 

Voo = lim — ^ 0. 

i^oo Sj 

Similar to case 1, Vqq must be of the form cx ■ in some half plane; here 
u is again a unit vector. On the other hand Proposition 14. 31 5) — which 
incidentally is preserved under the blow-up limit — implies that there 
are at least two vectors 7^ G dBi and two half lines {ae^ : a > 0} 
and {ae^ : a > 0}, such that ■ u > 0, Vvoo = Vfoo ■ e^, ■ u > 
and Vvoo = Vfoo ■ which is not possible for voo{x) = cx ■ u unless 
c = 0. □ 

Remark 6.5. Using methods in the forthcoming paper [TT], it is ac- 
tually possible to show that Case 1 in the proof of Theorem 16.41 cannot 
occur. More precisely, a degenerate point x° at which the set {vo = 0} 
contains a disk touching is not possible. 

Consequently the origin must in our example be a degenerate point. 
The fact that we can introduce many symmetry lines suggests that we 
can construct degenerate points with growth < {xl"^ in B1/2 for any 
integer m. 
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